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Abstract 

Finite difference solution of nonlinear partial differential 
equations requires discretisations and consequently grid errors 
are generated. These errors strongly affect stability and con- 
vergence properties of difference models. Previously such errors 
were analyzed by linearizing the difference equations for solu- 
tions. In this article properties of mappings of decadence [1,2] 
were used to analyze nonlinear instabilities. Such an analysis 
is directly affected by initial/boundary conditions. An algorithm 
has been developed, applied to nonlinear Burgers' equation [3,4] 
and verified computationally. A preliminary test shows that 
Navier-Stokes ' equation may be treated similarly. 
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1. The Objective , 

Let us consider a nonlinear partial differential equation 
8u/S t « L (u) (1.1) 

where L is a one-dimensional differential operator in x. Let the 
domain of integration be [a,b] X [0,«°). Equation (1.1) is subject 
to certain initial/boundary conditions and it is assumed that the 
problem is mathematically well-posed. 

An explicit finite difference analog of (1.1) is 

U n *» F(U n_1 ) (1.2) 

where, U n = (u” ••• u”) T £ D c R 1 , (R 1 = the real I-dimensional 

space), U? ■> iMx^t^) = the net function corresponding to u^ which 

is the true value of u at (x.,t ). 

i n 

An implicit finite difference analog of (1.1) is: 

G(U n ) = U r ‘ -1 . (1.3) 

Also, F: D c R 1 D and so is G. It is assumed that the trunca- 
tion errors are small and their effects are negligible. 

Grid error is defined by 

e 11 = u n - U n . (1.4) 

Stability is guaranteed iff V n J j|e n || < K, where K is positive and 
arbitrarily chosen. 

In this article an attempt will be made to see how one can 
obtain, V e 1 £ P. 1 

lim ||e n |[ = 0 (1.5) 

00 
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for given Ax (mesh size) and At (time step) . 

Obviously (1.5) guarantees stability. It also implies con- 
vergence for steady state solution. 


2 . Mathematical Preliminari es. 
Let, V n, z n t R 1 , and 


z n = A^"" 1 . (2.1) 

Clearly, lim z n = 0 iff 
n-*-* 

lim A A .... A « 0. (2.2) 

n n~ 1 1 

n->® 

Nov? (2.2) is true if there exists a particular norm such that 
V n > N 

II A J < a < 1 . (2.3) 

(These are discussed in details in 12].) Under these conditions 

(2.1) is said to describe a motion of decadence and A is called 

n 

a D-matrix. 

If instead of (2.1) the motion is given by 


A z 1 
n 


= z 


n- 1 


(2.4) 


it is a motion of decadance iff A 
if 


-1 

n 



< CL < 


1 


is 


a D-matrix which 


is true 
(2.5) 


for some particular norm and V n > K. 

It may be proved: 

Theorem: 1 If is a lower triangular matrix and 

p (A ) < a < 1, A is a D-matrix. (p(A ) = Spectral Radius of A .) 

n n n n 
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Theorem: 2 If A n is a tridiagonal matrix and (i) for i / j, 


a?J < la?. 


and (ii) ! a? ± - (a^ i _ 1 / ./aj^ # i _ 1 ) | >1, lajj 


A 1 is a D-matrix. The same is true if A n is a bidiagonal matrix 
n 

with nonnull elements on the main diagonal. 


3 . Analysis of Discretization Errors , 
Let us consider (1.2). Let 


_ # n~ 1 1 __ i .. n h _ n ™1 

F (u ) - F (U ) = A e 

n 


(3.1) 


Obviously, if a n . is an element of A , a? . = a?.(u n ,U n ) 
J ID n id iD 

the grid error equation for (1.2) is: 


Then 


e n = A e 11 " 1 . 
n 


(3.2) 


Hence, (1.5) is true if A is a D-matrix. 

n 

If we express, 


G(u n ) - G(U n ) = A e n 

n 


(3.3) 


then for (1.3), the equation (1.5) is true if A" is a D-matrix, 

n 

It may be seen that the effects of truncation error are total - 
ix. neglected in this discussion. Such effects were discussed in 


Thus, for an explicit finite difference equation, grid error 

effects are damped out if A in (3.2) is a D-matrix; and for an 

n 

implicit finite difference equation, the same is true if in 

(3.3) is such that A~ 1 exists and is a D-matrix, 

n 


4. Algorithm for Stability Analysis. 


It is well known that for any square matrix A (I X I) 
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(4.1) 


i n 

max a. . 

ij 13 

for certain natural 
(1.2), (1.5) is true 


I * |l\ll < I -max |ajj| 

norms. Thus, for an explicit equation like 
if 


I -max |a?.| < a < 1. (4.2) 

il 

If in case A is a lov?er triangular matrix, Theorem: 1 may be 
n 

applied. 

For an implicit equation of the form (1.3), if A is a tridi- 
agonal matrix, grid error effects may be studied by using Theorem: 
2. A general analysis for A^ (or A^ 1 ) to be a D-matrix may be 
found in [5]. 


5 . Application . 

Let us consider the inviscid Burgers' equation: 

u t + (1/2) (u 2 )^ =0. (5.1) 

Let the initial conditions be: 


u(x,C) = V if x < x (5.2) 

1 J 

=V 2 lfX> V 

V >V 2 

Let u be approximated by a tv;o-point forvzard difference for- 
mula and (u ) be approximated by a two-point, backward difference 
formula. Then the difference approximation of (5.1) is: 


n+l , n .2 i nv2 ( n, n 

u i = a(u ) ~ a(u ) + u i + T i 


If u* is replaced by U^* and x^ (the truncation error) is dropped, 
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then using e” = u” - 


U“, ve get: 


.n + 1 


= a (u n 


U n )e n + 
i-li-1 


{1 - a(uj + U”) }ej (5.3) 


where a = At/ (2Ax) . 

The linearized stability analysis requires: 

a(2V 1 ) < 1 (5.4) 

where V = max |u”|. This inequality implies restriction on time 
1 i f n 1 

step given by: 

At < Ax/V 1 . (5.5) 

In the present analysis (5.3) may be expressed as: 


e n41 = A e n (5.6) 

n 

where A^ is a bidiagonal matrix having diagonal elements a^ “ 

1 - a(u n + U?) and elements below the main diagonal as a? . , = 

ii 3 i, i -1 

a(u? , + u” ,) . Then by Theorem: 1, A is a D-matrix if 
l-l i+l J n 

max la” . 1 s a < 1, Vn>N. (5.7) 

1 1 1 1 

i 

If one chooses arbitrarily = 1.3, = 0.0, At = Ax = 0.1 

(and x T = x.), the linearized stability criterion (5.5) is violated, 
although (5.7) is satisfied. Computationally, instabilities were 
not found and the results given by fig. 1 seem to be quite cor- 
rect. 

Stability analysis of other explicit finite difference analogs 
may be treated similarly or by using the inequality (4.1). 

If both u t and (u 2 ) are approximated by two point backward 
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IIP I ''■! 


U 


difference formulas, we get an implicit finite difference analog of 
(5.1) and dropping the truncation error, the grid error equation 
becomes : 


• a K-i + + + + - e i _1 « 


(5.8) 


Here, is a diagonal dominant lov?er triangular matrix and 

la^j >1 V n > N. Hence, the numerical scheme is unconditionally 

stable by Theorem: 2. 

Let (5.1) be expressed as: 


u fc + uu x = 0. 


(5.9) 


If u fc is approximated by a two-point backward difference formula 
and u^ is approximated by a central difference formula, the error 
equation becomes: 


- au i e i-i + {i + + 


Here, is a tridiagonal matrix and considering the initial condi- 
tions (5,2), l a iil £ i* Hence, Theorem: 2 cannot be applied. Thus, 
stability criterion is not satisfied. (Linearized stability cri- 
terion is, however, unconditionally satisfied.) Actual computa- 
tions showed instabilities. Now if we change the initial boundary 
conditions as: u(x,0) - x, u(0,t) = 0, u(l,t) = 1/(1 + t) , 

u. _ > u. „ Vi and la 1 ?. I > 1 with diagonal dominance, the im- 
plicit scheme should now be unconditionally stable. ziebarth (6] 
verified it computationally. 


6 . A Remark on Navier- St okes T Equation . 

Let us consider Navier-Stokes’ equation in the vcrticity-streaiu 
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function form as: 


C t + ^y ' ^y^x = vv2c (6,1) 

V 2 i p = -Z (6.2) 

where ? = vorticity and ip = stream function. This coupled system 
is subject to some specified initial-boundary conditions. If we 
analyze the grid errors for implicit schemes we get two equations 
of the form 

* n e n + e n f n = e n ~ 1 (6.3) 

A f n = e n (6.4) 

n 

where e n = grid error for £ and f n = grid error for ip [7]. These 
equations may be expressed as 

A n e n = e n_1 . (6.5) 

It appears that if sharp discontinuities are present neither in the 
flow field nor on the boundary, conditions of Theorem: 2 will be 
satisfied. Therefore, the implicit scheme will be stable. 

7 . Conclusion . 

If the sequence of matrices {A n } be such that V n, 

|1 A |j < a < 1, ||e n || will form a monotone decreasing sequence, where- 
as if V n > N, [Ja^II < a < 1, ||e n || may show some oscillations before 
it is damped out* In both cases, however, as n ■+ » || e n j| 0 • 

For the linearized grid error analysis, A n = A V n and if A is 
a convergent matrix stability is obtained. Thus, linearized grid- 
error theory is a particular case of the analysis presented here. 


Since elements of A are functions of u n and U n , initial- 

n 

boundary conditions affect the properties of A n . 

In order to check that A (or A -1 ) is a D-matrix, some inf or- 

n n 

mation regarding the nature of the solution must be known a priori. 
This may be done mathematically or experimentally or both. 
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Figure 1 


Ax = At = 0.1 

Shock speed = 
(V 1 + V 2 )/2 

o 


“ Initial solution 
O Solution at 0.4s 
A Solution at 0.8s 
□ Solution at 1.2s 
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* Explicit finite difference solution of equation (5.1). 
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